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Abstract 

In the ZEUS experiment at HERA, photoproduction processes have been studied for photon-proton 
centre-of-mass energies in the range 100 < W lp < 300 GeV and jet transverse energies extending to Erf? ~ 
70 GeV. The data contribute to our understanding of QCD dynamics, and also provide new constraints on 
the photon's parton density. 
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1 Introduction 



There has been tremendous progress in the understanding of photon induced reactions since the start of HERA 
operation. Up to the end of 1994, HERA had delivered a total of 7 pb _1 of ep data. The analysis of this data 
was fruitful. ZEUS was able to clearly establish evidence for both the direct and the resolved classes of process. 
These are illustrated at leading order (LO) in Figure's 1(a) and (b). 




Figure 1: Examples of leading order diagrams for (a) direct photoproduction, (b) resolved photoproduction and 
(c) direct prompt photon photoproduction. 



An observable quantity was defined which extended the distinction between direct and resolved processes to 
all orders of perturbation theory, x® BS — JDjets E^e -71 ' /2E 7 where the sum runs over the two highest trans- 
verse energy jets in the final state. Direct processes contribute primarily at high a:^ B and resolved processes 
contribute at low x^ BS . However the comparison with theory has been hampered by the need to introduce 
models such as secondary interactions between the photon and proton constituents in order to understand the 
resolved processes at ~ 6 GeV. 

By the end of 1997 running HERA had delivered an order of magnitude more data, - 70pb~ 1 . This large 
data-set allows the exploration of photoproduction up to much higher jet transverse energies where a strong 
confrontation with the predictions of perturbative QCD (pQCD) may be made. This is illustrated in Figure ||| 0. 
Here the x® BS distribution is shown for photoproduction events where at least two jets are produced, one with 

> 14 GeV and the other with > 11 GeV where the jet pseudorapidities satisfy -1 < ?7 jot < 2. The 
results are compared with the Monte Carlo predictions from HERWIG which implements the leading order 
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Figure 2: Distribution of x° BS . The black dots show the uncorrected data. The histogram shows the HERWIG 
prediction with the contribution due to LO direct processes shaded. Multiple interactions have not been included 
in the simulation. 

matrix elements together with a parton shower. A good description of the data is provided by HERWIG across 
the entire x^ BS range, without having to resort to additional parameters. This facilitates the interpretation of 
the data within pQCD. 

This high accumulated luminosity also gives access to such infrequent processes as prompt photon production, 
illustrated in Figure 0(c), and multijet production. In this talk results are presented on four topics; single and 
dijet inclusive cross sections |l| [jj, prompt photon production jj], the evolution of the photon's structure with 
virtuality || and two and three-jet angular distributions in high mass events 0, [?]]. 



2 Single and Dijet Inclusive Cross sections 



In hadron-hadron collisions, jets of final-state hadrons are commonly determined by cone algorithms which 
seek to maximize the transverse energy produced inside a cone of fixed radius in pseudorapidity and azimuth. 
However the confrontation of data and theory using such algorithms is compromised by the need to introduce 
an ad-hoc parameter -Rsep into the theory, in order to simulate merging effects S. A clustering algorithm 
which combines particles with small relative transverse momenta, kr, into jets has been developed for use in 
hadron-hadron collisions ||, . With this algorithm it is not necessary to introduce additional parameters in 
order to compare data and theory. In the results presented here the kx clustering algorithm has been used. 

Figure || presents the inclusive jet cross section, da/dEj!, for jets with —0.75 < Tf ct < 2.5 and for centre- 
of-mass (CM) energies in the range 134 < W lp < 277 GeV M. The cross section falls by 4 orders of magnitude 
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Figure 3: Left: Inclusive jet production cross section, da/dE^ T . Right: Fractional difference between the 
measured cross section and the pQCD prediction using the GRV HO parametrization of the photon parton 
densities. The dots show the data. The inner error bars show the statistical error. The outer error bars show 
the quadratic sum of the statistical and systematic uncertainties with the exception of the absolute energy scale 
uncertainty which is shown separately as a shaded band. 

over the E^ range, extending to ~ 70 GeV. This behaviour is impressively reproduced in the theory. The three 
curves show 0{aa 2 s ) pQCD calculations for different parametrizations of the parton densities of the photon. 
Within the current measurement uncertainties all are consistent with the data. This is more clear in the plot on 
the right which shows the relative difference with respect to the calculation using the GRV HO parametrization. 

The dijet inclusive photoproduction cross section as a function of the pseudorapidity of one of the jets, 
da/dr]2, is shown in Figure ||, in three bins of the pseudorapidity of the other jet, 771 Here the events have 
been required to have at least one jet with E^ > 14 GeV and at least one other with E 3 £ > 11 GeV. The 
photon-proton CM energies have been limited to the range 212 < W lp < 277 GeV, in order to enhance sensitivity 
to the photon's parton densities. Indeed, by comparison with the next-to-leading order (NLO) calculations we 
find that one parton density set (GRV) is favoured over another (GS96). It is worth emphasizing that GS96 is 
a modern set, consistent with the world data on F% . The ability of the data to disfavour this, indicates that 
ZEUS is contributing photoproduction information in a heretofore unexplored regime. 




3 Prompt Photon Production 

A clean laboratory for testing QCD predictions is provided by prompt photon production, Figure 0(c). An 
inclusive prompt photon measurement, tagging only the final-state hard photon, may be directly compared 
with pQCD calculations, with neither concern for jet definition matching between theory and experiment, nor 
for hadronization uncertainties. 

In Figure |B| the inclusive prompt photon cross section with respect to the photon's pseudorapidity, da jdrp , 
is shown. The events lie in the range 120 < W lp < 270 GeV and have at least one photon of transverse energy 
5 < El < 10 GeV. The cross section is compared with predictions for two different parametrizations of the 
photon's structure, GRV and GS. Again we find that there is a preference for the GRV parton densities over 
the GS parton densities. 
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Figure 4: Dijet inclusive cross section, da/drf^*. Error bars are as described for Figure^. NLO calculations for 
two parametrizations of the photon parton distributions are shown by the solid and dashed curves. 
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Figure 5: Inclusive prompt photon photoproduction cross section, da/drp . The errors shown are statistical 
only. The lines show QCD NLO predictions for two different parametrizations of the photon parton densities. 
A systematic uncertainty of ±15 % should be allowed. 



4 Virtual Photon Structure 

In addition to asking about the parton density of the real photon, one can ask how that structure will evolve 
with the negative invariant mass squared of the photon, or the photon's virtuality, Q 2 . One expects that as Q 2 
increases, less time is available for the photon to develop a complicated hadronic structure. 

In Figure ^, the uncorrected distribution of x° BS is shown for events with at least two jets of > 6.5 GeV 
and —1.125 < rf et < 1.875, in three different ranges of Q 2 . These distributions indicate that the frequency 
of resolved processes (which populate low x® BS ) decreases with respect to the frequency of direct processes, 
as Q 2 increases. To quantify this observation, the ratio of the dijet inclusive cross section for x® BS < 0.75 to 
the dijet inclusive cross section for i° BS > 0.75 has been measured as a function of Q 2 . This is also shown in 
Figure |^. This measurement confirms the expectation that resolved photoproduction processes are supressed as 
the photon's virtuality increases. 
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Figure 6: Uncorrected x™ eas (a;° BS ) distributions in three photon virtuality ranges where the error bars show the 
statistical errors only. The histogram shows the HERWIG prediction with the LO direct contribution shaded. 
Also, the ratio of the resolved to the direct dijet cross sections versus Q 2 where the error bars are as described 
for Figure ||. 



5 High Mass Dijet and Three-jet Distributions 

As well as probing the partons within the photon, the ZEUS photoproduction data can be used to learn 
something about the dynamics of parton-parton scattering. The definition of the scattering angles in two and 
three-jet events are illustrated in Figure |?]. 




Figure 7: Centre-of-mass frame diagrams of two and three-body scattering angles. The beam direction is 
indicated by Pg- 



In Figure g the dijet invariant mass distribution is shown for events in which the scattering angle in the 
dijet rest frame, 6*, satisfies | cos 9* \ < 0.8 M. Resonances decaying to two jets would be expected to show an 
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Figure 8: Left: Differential dijet cross section, da/dM 3J . Right: Differential dijet cross section, da/d\ cos 9* \. 
The error bars are as described for Figure |^. NLO QCD calculations for two parametrizations of the photon 
parton distributions are shown by the dashed and dot-dashed lines. 



excess above the pQCD predictions in a limited range of M JJ . No evidence for this is seen. 

The scattering angle in the dijet CM frame in photoproduction has been shown to be sensitive to the spin 
of the exchanged parton |ll| . The cos 9* distribution is shown for the events with Mj j > 47 GeV on the right 
in Figure ||. The NLO pQCD calculations are in good agreement with the data confirming that the dynamics 
of two-jet production are well understood. 

ZEUS has also measured the high-mass three-jet cross section, da/dM^j, as shown in Figure || 0. The 
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Figure 9: The three-jet cross section da/dM 3J The error bars are as described for Figure [|. 0(aa 2 ) pQCD 
calculations by two groups of authors are shown by thick solid and dot-dashed lines. The thin solid and dashed 
histograms show the predictions from PYTHIA and HERWIG. 

0(aa 2 s ) pQCD calculations from two groups of authors (Harris & Owens, Klasen & Kramer) provide a good 
description of the data, even though they are leading order for this process. Monte Carlo models also generate 
three-jet events through the parton shower mechanism and both PYTHIA jl^] and HERWIG reproduce the 
shape of the M 3J distribution. 

For three-jet events there are two relevant scattering angles as illustrated in Figure 0. The distributions of 
cos #3 and ipa are shown in Figure O. The cos #3 distribution resembles that of cos 8* and exhibits forward and 
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Figure 10: Distribution of the angles cos 6*3 and ip3- The dotted curve shows the distribution for a constant 
matrix element. Other details are as described for Figure ||. 

backward peaks. It is well described in both (D(aa 2 s ) pQCD calculations and parton shower models. The ip3 
distribution is peaked near and n indicating that the three-jet plane tends to lie near the plane containing 
the highest energy jet and the beam. This is particularly evident if one considers the ^3 distribution for three 
partons uniformly distributed in the available phase space. The phase space near ^3 = and n has been 
depleted by the cuts and by the jet-finding algorithm. The pQCD calculations describe perfectly the ^3 
distribution. It is remarkable that the parton shower models PYTHIA and HERWIG are also able to reproduce 
the ip3 distribution. 



Within the parton-shower model it is possible to determine the contribution to three-jet production from 
initial-state radiation (ISR) and final-state radiation (FSR). In Figure |ll| the cross section da /dips is shown, 
along with the area-normalized ip3 distribution. Here the PYTHIA prediction is shown for FSR only, ISR only, 
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Figure 11: (a) The cross section da /dips and (b) the normalized distribution of ips- The error bars are as 
described for Figure^. The solid line shows the default PYTHIA prediction, the dashed line shows the prediction 
with FSR switched off and the dot-dashed line shows the distribution with ISR switched off. 



and default PYTHIA which includes the interference of these two. Clearly, ISR is predominantly responsible 
for the three-jet production. 

The QCD phenomenon of colour coherence may be expected to favour topologies in which the softest jet 
is radiated close to the plane containing the highest energy jet and the beam. Figure 12 shows da /dips and 
1/ 'ada j 'dips compared with the expectations from HERWIG, default PYTHIA, and PYTHIA with coherence 
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Figure 12: (a) The cross section da/dips and (b) the normalized distribution of ips- The error bars are as 
described for Figure ||. The PYTHIA and HERWIG predictions are shown by the solid and dashed lines and 
the dashed-dotted line shows the PYTHIA prediction with colour coherence switched off. 



switched off. Coherence suppresses large angle emissions leading to a depletion of the ips distribution near 
ip3 = tt/2. The incoherent prediction is relatively flat. Coherence in the parton shower is required by the data. 
However the current measurement precision is insufficient to discriminate between the two models, HERWIG 
and PYTHIA. 



6 Summary 

ZEUS has measured hard photoproduction in a kinematic regime where the photon's structure has not yet 
been well-constrained. In both inclusive prompt-7 production and in dijet production the GRV photon parton 
densities are favoured over the GS96 parton densities. The data reveal fundamental phenomena of the theory 



of QCD. A suppression of the resolved component of the photon with increasing photon virtuality has been 
observed and the effect of QCD dynamics on the angular distribution of two- and three-jet events has been 
studied. New photoproduction territory will shortly become accessible to the HERA experiments as we move 
into a regime of high luminosity HERA running. 
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